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Abstract 

Considering the fact that some excited states of the heavy quarkonia (charmonium and bot- 
tomonium) still missing in experimental observations and potential applications of the relevant 
wave functions of the bound states, we re-analyze the spectrum and the relevant wave functions of 
the heavy quarkonia within the framework of Bethe-Salpeter (B.S.) equation with a proper QCD- 
inspired kernel. Such a kernel for the heavy quarkonia, relating to potential of non-relativistic quark 
model, is instantaneous, so we call the corresponding B.S. equation as BS-In equation throughout 
the paper. Particularly, a new way to solve the B.S. equation, which is different from the tradi- 
tional ones, is proposed here, and with it not only the known spectrum for the heavy quarkonia is 
re-generated, but also an important issue is brought in, i.e., the obtained solutions of the equation 
'automatically' include the 'fine', 'hyperfine' splittings and the wave function mixture, such as 
S — D wave mixing in J PC = 1 states, P — F wave mixing in J PC = 2 ++ states for charmonium 
and bottomonium etc. It is pointed out that the best place to test the wave mixture probably is 
at Z- factory (e + e~ collider running at Z-boson pole with extremely high luminosity). 
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I. INTRODUCTION 



Spectroscopy, including the spectrum and the corresponding wave functions, is a very in- 
teresting topic for heavy quarkonia in particle physics. The spectrum and the corresponding 
wave functions for the binding systems can be tested experimentally and via study of the 
spectroscopy one may have insight of the heavy quarkonia and understand QCD , which is 
the nature of the binding, further as well. In the literature, there are various approaches to 
the spectroscopy of the heavy quarkonia: charmonium and bottomonium 1KZ|, and to solve 
the Bethe-Salpeter (B.S.) equation is one of them 4|-[7|]. Since recently we have realized a 
new method to solve the B.S. equation for the heavy quarkonia, so in this paper, we would 
like to try the method i.e. to apply this method to re-analyzing the spectroscopy of the 
heavy quarkoia: charmonium and bottomonium under B.S. equation approach. 

First of all, how to determine the B.S. kernel is crucial for B.S. equation approach to 
a bound state problem. It is known that if one adopts the QCD-inspired Bethe-Salpeter 
(B.S.) equation [8j] approach to the problems of hadronic bound states, then the relevant B.S. 
kernel for a double heavy quark- ant iquark system, such as charmonium and bottomonium, 
is instantaneous approximately i.e. the B.S. equation is essentially an instantaneous one (a 
BS-In equation). It is also known that BS-In equation can further precisely relate to the 
Schodinger equation in potential model (PM) by means of the Salpeter approximate method 
9]. Therefore, one may use the relation to the potential model (PM) and help oneself to 
determine the kernel of BS-In equation precisely. Whereas starting with the BS-In equation 
whose kernel is fixed in terms of QCD consideration and the relation to PM, one can extend 
some relativistic nature of the problem more than what PM can consider, and stand on 
more solid theoretical ground for the B.S. equation approach, hence, we start the study of 
spectroscopy for heavy quarkia with such a BS-In equation. Moreover, we apply the new 
realized method to solving the BS-In equation. People later on will see an important issue 
from the new method is that besides the 'fine and hyperfine' splitting being involved, the 
wave mixtures in the wave functions, such as S — D wave mixing in J = 1 (J: total 
angle momentum; P: parity; C: charge parity) states and P — F wave mixing in 
states etc, are determined precisely, although the mixtures, in fact, are rooted in the kernel 
under the present framework of BS-In equation. 

The new proposed method can be outlined as that, firstly we analyze the bound states 
according to their total angular momentum J, parity P and charge conjugation C, such as 
the states (T + (%), l~"( 3 Si or 3 D 1 ), ++ ( 3 P ), l ++ ( 3 ^i), 2++( 3 P 2 or 3 F 2 ), and etc; 
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secondly we write down the most general formulation for the B.S. wave functions respectively, 
and then input the formulated wave functions into the BS-In equation and turn the equation 
into a set of proper coupled equations for the components which appear in the formulation; 
finally we solve the coupled equation numerically, and obtain the mass spectra and wave 
functions for (cc) and (bb) binding systems. For convenience, we call the coupled equations 
as BS-CoEqs later on. 

This paper is organized as following, additional to the Introduction section I, in section II 
we introduce the relativistic Bethe-Salpeter equation and BS-In equation. In section III we 
start with the generalized formulation for relativistic wave functions with definite quantum 
numbers to derive the relevant BS-CoEqs for low total angle momentum states individually. 
Finally, we show the numerical results although we do not present the detail to solve the 
equation numerically, and we also explain and briefly discuss the obtained solutions of the 
BS-CoEqs in section IV. 



II. INSTANTANEOUS BETHE-SALPETER EQUATION 

Firstly let us outline the reduction of the B.S. equation which is similar to the way of 
Salpeter [9[ if B.S. kernel is instantaneous, and introduce necessary notations. The readers, 
who are interested in the details, can also find them in Ref. 1CH12|. 

The Bethe-Salpeter (B.S.) equation for mesons is read as: 

{A - m 1 ) X {q){^2 + m 2 ) = ij 7^4v(P, K q) X (k) , (1) 

where x{q) is the B.S. wave function, V(P, k, q) is the interaction kernel between the quark 
and antiquark, and pi,p 2 are the momenta of the quark 1 and anti-quark 2. Quark mass 
is mi, antiquark mass is m 2 , and here we consider heavy quarkonia: charmonium and 
bottomonium so we have nil = m 2 . The total momentum P and the relative momentum q 
are defined as: 

mi 

Pi = aiP + q, ai 



p 2 = a 2 P - q, a 2 



mi + ^2 
m 2 



nil + m 2 

and «i = «2 = \ for charmonium and bottomonium. 

We divide the relative momentum q into two parts, q\\ and q±, 
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gf = (P ■ q/M 2 )P» , = tt ~ 4 
Correspondingly, we may have two Lorentz invariant variables: 



When P= 0, they turn to the usual component q and |g| respectively 
If the kernel V(P, k, q) takes the simple form: 

V(P,k,q)^V(k ± ,q ± ) 

namely the B.S. equation is 'instantaneous', for convenience, we would like to introduce the 
notations (f p (q±) and r](q^_) so the 'instantaneous (three dimensional) objects' will accord- 
ingly read as follows: 

<pA4D = i[^rx(4!,4D, (2) 

l{<£) = J -^V(k ± ,q ± )<p p (ktt . (3) 
The B.S. equation now is rewritten as: 

x(q\\,Q±) = S 1 (pi)ij(q±)S 2 (p2) ■ (4) 
Generally the propagators of the two constituents can be decomposed as: 

SM = a+ + Kp(q±) ; (5) 

' J{i)q. P + a iM — 0Ji + ie J{i)q P + a^M + Ui — ie ' 



with 

1 

2a;,- 



u i = Jm 2 i +q 2 Kf P { 



P 

—Ui ±J(i)(mi + 4 ± ) 



(6) 



where i = 1, 2 for quark and anti-quark respectively, and J(i) = (— and Af P (q±) satisfy 
the relations: 

A£>(?j_) + Kp(Q±) = ^ , A l ± P (g ± )^Af P (g^) = A± (q ± ) , A± (^A^fa.) = . (7) 

Hence sometimes Af P (q±) are called as 'project operators', although they need to be sand- 
wiched with the operator jj when 'projecting' as Eq©. 

Introducing the notations ip F z± (q±) to note the projected wave functions as: 

^(q ± ) = AtpM^AlD^pM , (8) 
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and we indeed have 



With contour integration over q p on both sides of Eq.(j4]), we obtain: 

Ajp(<a)r/(g ± )A+ P (g ± ) K P {q±)ri{q±)h.2p{q±) 



M9J 



(M - 2wi) (M + 2wi) 

and the equation becomes four independent equations: 

(M - 2u 1 )i f + + (q 1 _) = At P (q ± Mq±)Atp(q±) , 
(M + 2wi)^;-(gj.) = -Aip{q±)v{q±)A2p{q±) , 

vt~(9d = v~ + (q±) = • (9) 

where we have u\ = U2 for the equal mass system. In fact the four equations is of an 



'eigenvalue problem' about the eigenvalue M. Note that in the Ref. [9| the way for solving 
the BS-In equation is not exactly equivalent to the four equations Eq.flH]). Details about 



examining the equivalence may be found in Ref. 1(| . Alternately here we exactly start with 
the four equations to solve the BS-In equation. 



The normalization condition for B.S. wave function is: 



q T dq T 
2tt 2 



Tr 



M M AT M 



2P ■ (10) 



Now let us return to the problem for the heavy quarkonia (cc) and (66). To fix the kernel 
for the heavy quark and heavy anti-quark, on one hand, we should let the kernel being 
QCD-inspired and on the other hand, we should relate the kernel to the Cornell potential 
accordingly. Thus the kernel in space-time looks like as a linear scalar interaction (the 
confinement one in QCD nonperturbative nature) plus a vector interaction (single gluon 
exchange in Coulomb gauge): 

4 o 

V(r) = V a (r) + 7o ® 7°K(r) = Ar + V - % <g> 7%— , (11) 

6 r 

where A is the string constant, a s (r) is the running coupling constant of QCD. Usually, in 
order to fit the data of heavy quarkonia, a constant Vo is often added to the scalar confining 
potential and takes different values for the bound states with different quantum numbers 
respectively 1 . 



1 One will see later on in this paper that the value of Vo is determined by fitting the data for the ground 
states with the corresponding quantum numbers. 
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To avoid the infrared divergence in the Coulomb-like one and to correspond the fact that 
the confined linear interaction should be also suppressed at large distance phenomenologi- 
cally, so it will be better to re-formulate the kernel as follows: 

V s (r) = ±(l-e- a ' r ) + V , 

V v (r) = -lfe~ ar . (12) 

To decrease the parameters which are needed to fix by fitting data, we assume a' = a 
approximately 2 . It is easy to show that when ar < 1, the potential approximately becomes 
linear. Now the B.S. kernel in momentum space and in the rest frame of the bound state is 
read as: 



V(q) = V s (q)+ % ® 1 °V v (q) 
V.CS) = -fi + V Q )P@ X 1 



vjt)= 2 



37r2 ( f + a 2 ) 



/-N 1271 1 
a.(g) = - — (13) 

66 — Zl\ f i i q \ 

1 log a + 

QCD 

where Nf = 3 for (cc) system, Nf = 4 for (66) system; the constants A, a, a, Vq and Aqcd 
are the parameters which characterize the kernel (potential). 



III. GENERAL FORMULATION FOR THE B.S. WAVE FUNCTIONS AND THE 
COUPLED EQUATIONS 

In fact in this section, we show the new realized method to solve a BS-In equation, but 
specifically apply to the concerned heavy quarkonium problem. 

Firstly, according to the total angle momentum (J), parity (P) and charge conjugation 
(C) of the concerned bound state, we write down the most general formulation for each of 
the relativistic B.S. wave functions, and then we put it into Eq.(j9]) to derive out the coupled 
equation for the components appearing in the formulation, BS-CoEq. In the below subsec- 
tions, we do the derivation for the low-laying states: J PC = 0~ + , 1 , l + ~, ++ , 1 ++ , 2 ++ , etc 
in turn precisely. 



2 In fact, at final step (numerical solving BS-CoEq) we find that the results are not very sensitive to the 
assumption when a and a' vary in reasonable region. 
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A. J PC = + , 1 + , ++ and 1 ++ wave functions and BS-CoEqs for relevant com- 
ponents 

Since the bound states with the quantum numbers J PC = _+ , 1 H , ++ and 1 ++ are 
similar, so in this subsection we derive the equations for them in turn. 

I. The bound states with quantum numbers J PC = 0~ + , which in non-relativistic frame- 
work are 1 Sq states mainly. 

The general formulation of the In-BS wave function Eq.(12T) for the states J PC = 0~ + is 

HQ 



Pfi(q T ) + Mf 2 (q 7 



M 



1.. 



5 ' 



(14) 



M2 



where M is the mass of the bound state (the corresponding meson) and = q^ 
is the four dimensional vector. In the center mass system q± = (0, <?), q T = | Q \. Later on 
we abbreviate q T as q if it does not make any confusion. 

Now let us derive the coupled equations (BS-CoEqs) from Eq.flH]). To put Eq.(fl4l) into 
the last two equations of Eq. ([9]) 

^-M = ^-M = ° 

and by taking various traces for 7-matrices on both sides of the equations, we obtain the 
independent constraints on the components for the wave function: 

fi(q)M 



(15) 



so we can apply the obtained constraints Eq. ffl5|) to Eq. ffl4|) and rewrite the relativistic wave 
function of state 0~ + as: 



mi 



[l + —)rfi{q) + Mf 2 (q) 
mi 
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(16) 



From the above formulation of the wave function one can see clearly that besides the 'great 
component', which is proportional to either M75 or ^75, there is also a 'small component', 
which is proportional to ^^75, linear in q± (P-wave nature) and suppressed by 

Put the wave function Eq.f lTB]) into the first two equations of Eq.()H]) and by taking various 
traces for 7-matrices to both sides of the equations, we obtain the independent coupled 
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integral equations (BS-CoEqs): 



(M - 2u>i) 



A(ff) + /a(9) 



mi 



d 3 k 1 
(2tt) 3 u;? 



x{{V s -V v ) fiityml + MtymxUi - (V s + V v )fi(k)(q ■ k) 



(M + 2ut 



/i(?)-/ 2 (g)— 



d 3 k l 

(2tt) 3 u;? 



x (K-K) fiifyml-fcifym^i - (V s + V v )fi(k)(q ■ k) 



(17) 



here = Wi = ynii + q^- Now we are prepare ready to solve the BS-CoEqs Eq. ffTTl) . 
as an eigenvalue problem, for fi and f 2 numerically, specially in center mass system, and 
we may obtain the required results (the spectrum for J PC = 0~ + states and the B.S. wave 
functions accordingly) finally. 

Now accordingly the normalization condition is read as 



d 3 q 



4A(g)/ 2 (g)M 2 



OJi mi 



+ — + 



2M . 



(18) 



(27r) 3 [mi oji ojitrii 

II. The bound states with quantum numbers J PC = l + ~ which in non-relativistic frame- 
work are 1 Pi states mainly: 

As that for the states J p = 1 H , the general form of the In-BS wave function can be 



written as 



12 



16|: 



PA+- 



From the equations 



(q) = q±-e X ± h{q) + f 2 (q)^ + fM^t 



75- 



(19) 



(20) 



a constraint on the components of the wave function 

f*(q)M 



h{q) 



nil 



is obtained. 

With the constraint, the wave function now turns into: 



p,i+- 



(?) = q±-t } _ 



A(q)+f2(q)£-h(q)^r 
M mi-M 



7s 



q± -e_ 



/i(g) + /a(ff)(l + 



A) I. 

mi' M 



75 



(21) 
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here the factor (q± ■ e*) indicates the wave function is of P-wave nature mainly; whereas in 



Eq.( l2~Tj) the 'small component' term (q± ■ ei)^-^7s contains high order wave. 
Now the normalization condition for the 1 P\ wave function is read as: 

J (2tt) 3 3m 1 v ; 

In terms of the same derivation as that for J PC = h states, we obtain the coupled 
equations (BS-CoEqs) for the components f\ and f 2 : 

d 3 k (k ■ q) 



(M - 2cji 



/i(?) + /a(?) 



mi 



(27r) 3 Uirriiq 2 



x {(V s + V v )f 2 (k)(k-q) - (V s - V v ) [fi(k)u imi + f 2 (k)ml]} , 
(M + ^)f/ 1 ( ff )-/ a ( ff )^l=-/ 7 ^iH 

x {-(v s + v v )f 2 (k)(k ■ q) - (v s - v v ) - / 2 (>K] } 



(23) 



Therefore, we are prepare ready to solve the coupled equations BS-CoEqs, as an eigenvalue 
problem, for f\ and f 2 numerically, specially in center mass system, and we may obtain the 
required results (the spectrum for J PC = l + ~ states and the B.S. wave function accordingly) 
finally. 

III. The bound states J PC = ++ and 1 ++ which in non-relativistic framework essentially 
are 3 Pq and 3 P\ states respectively: 

Since the bound states J PC = ++ and 1 ++ are very similar, thus here we treat them 
simultaneously. For the states J p = ++ , the general form of the In-BS wave functions can 
be written as [ill [l6]: 



¥?++(<?) = h{q)4± + h{q) 



— + Mq)M • 



With the equations 



C(?) = C(!) = 



(24) 



(25) 



we obtain the constraints: 



Then the wave function: 



1M 



fi(q)q 2 
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P.0+ + 



(?) = /i(?)^L + /a(?) 

4x 



ix. P fl(q)q 



M 



mi 



mi M 



(26) 



Here the factor tf± = (q- 7) in CMS is contained in the wave functions, that means the the 
wave function is of P-wave nature. Whereas the term which contains q 2 = in 
Eq. (l26l) is suppressed by the factor --j. 

In terms of the same way as that for J PC = 0~ + states, with the first two equations of 
Eq.Q we obtain the coupled equations for the J PC = ++ states: 

d 3 k 1 



(M - 2ui) 



(M + 2wi 



fi{q)+h{q) — 



x V, + V v 



(2tt) 3 w 1 2 

+ ( mi )(V s - V v ) [fi(k) mi + f 2 (k)u 1 ] (k ■ q)} 



k 



fM-h{q) 



mi 



d 3 k 1 r r 



(mi)(K - V v ) [h(k) mi - h(k)ui] (k-q)} 



The normalization condition for the wave function is read: 

d 3 q 4/ 1 / 2 w 1 g 2 



/ 



(27T)* 



M . 



mi 



(27) 



(28) 



Whereas for the J PC = 1 ++ states, the general form for the wave function can be written 

as 



<P Pil++ (q) = ie^P v ql^[h{q)M^ + / 2 (g)/V 
+th(q)e wa5 q± p P,l8l 5 /M] /M 2 . 



From the he equations 



<; ++ (?) = < ++ (s) = o, 



(29) 



(30) 



we obtain the constraints on the components of the wave function: 

/ 2 (g)M 



/s(?) 



nil 



Then we have: 



<P P>1++ (q) = ie^ aP P v qt^[h{q)M^ + f 2 (q)(PY + ie^ S q^lsl,/ mi)]/ M 2 . (31) 

Here the front factor e pva pP h 'qf_e^ ', being linear in q, means the wave functions are of P-wave 
nature. 

In terms of the same way as that for J PC = 0~ + states, with the first two equations of 
Eq.fjH]) we obtain the coupled equations (BS-CoEqs) as follows: 

d 3 k 1 



(M - 2u x ) 



fx(q) + f2(q) — 
nil 



(2?r) 3 2uj 1 m 1 k 2 q 2 
10 



X {-(Vs + V v )f 2 (k) [k 2 q 2 + (k-q) 



(M + 2wi 



2m x {V s - V v ) [fx{k)uJx + f 2 (k)nn} k\k ■ q}} 
r db 1 r 

-/(W^¥ x { (l, ' + K)AW 

The normalization condition for the J PC = 1 ++ wave function is read: 



fM-h{q)— 

mi. 



d 3 q 8hf 2 u iq 2 



M. 



(2tt) 3 3mi 

Now we are ready to solve the coupled equations Eqs. (j27l 152]) numerically. 



k 2 q 2 + {k-q) 2 

(32) 

(33) 



B. J PC = 1- , 2 ++ wave functions for In-BS equation and BS-CoEqs for relevant 
components 

As for the states J PC = 1 ,2 ++ , they are quite different from the states in the above 
subsection, because there is S — D wave mixing in the J PC = 1 states and there is P — F 
wave mixing in the J PC = 2 ++ states. 

I. The bound states J PC = 1 which in non-relativistic framework are 3 Si and/or 3 Di 
states mainly: 

First of all, we write down the general formulation for the wave functions of In-BS equa- 
tion with quantum numbers J p = Y 
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151: 



(?) = q± ■ <d 



+A?h{q) + - Pq± ■ e X ± )f 2 (q) 



(34) 



where the is the polarization vector of the vector meson. From the last two equations of 
Eq.® 

(?) = v> x ;~l (q) = o , (35) 

we obtain the independent constraints on the components of the wave functions: 



A(?) 



-q 2 h{q) + M 2 / 5 (g) 



Mq) 



h{q)M 



Mni\ nil 
Then with the constraints, there are only four independent components /3(g), /4(g), /5(g) 
and fe(q) left in the Eq. (l34l) . Namely 



Vi-(gj-) = ?±-ei(- 



+ ^)Hq)+q±-ei^U(q) 



■Mmi M' 

+(mA + n ■ 4-)Mi) + [Ar + - ( - lAM ]fM, (36) 



mi 
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and from the formulation it is easy to realize that the 'great components' in the wave 
function, which are proportional to /s or / 6 and or (/[_ JP) are of S'-wave nature, whereas 
the components in the wave function, which are proportional to f 3 or f± and (q± ■ e\)4±_ 
or (q± ■ ei)^j_ JP (double q± being contained) are of D-wecve nature (a tensor about q±). 
Therefore, no matter what are the other 'small terms', the wave functions Eq. (l36|) involve 
S — D wave mixing properly. 

To put Eq. (l36|) into the first two equations of Eq.(j9]) and take various traces on both 
sides of the equations, we obtain four coupled integral equations for the four independent 
components fs, /s and /e (BS-CoEqs): 



(M - 2w a ) h{q) 



M 2 



h(q)) + [U(q)j^ + fe(q) 



mi 

LOi 



d 3 k 2 



(V s + V v )lf 3 (k)—-f 5 (k)\(k-q) 



~(Vs - V v ) 



ml h{k) 



(k-q)' 
M 2 q 2 



(k-q)' 
M 2 q 2 



+ U{k) 



(37) 



(M + 2ui)Uf 3 (q)j^-h(q) 



nil 



(V - V v ) 



m{ f 3 (k) 



d 3 k 2 
JZnfu 2 K 

(k-q) 2 



(V a + V) 



)-[fMj f2+ Mq))^ 

' fm^ 2 - h{k))\(k ■ q) 



M 2 q 2 



f 5 (k)\ -miui \ f 4 (k) { ±^ + f 6 (k) 
(M - 2u x ) { ( h(q) + h(q)—) TT, ~ 3 ( h(q) - /e(ff)^) - fe(q) 



d 3 k 1 



+(V S - V v ) 



(V + V 
k 2 



uoi) M 2 
2ui 



I 2 \ 



k 2 



< h(k) 



h(k) - h(k)— + h(k) 

nil Jm 



k 2 



niiUJi 
(k-q) 



M 2 



3f 5 (k)\+ mi ui[U(k)— + 3f 6 (k) 



M 2 



U(k) 



(M + 2ui) 



h(q)-U(q) 
d 3 k 1 



mi 



(2tt) 3 u; 2 



U)i 

(V + V v ] 



(k-q)' 
M 2 

l_ 

M 2 

2ui 



->2 



fc(k)q 



3 (Mq) + Mq) 



mi 



h(k)-h(k)— + h(k) 
mi M z 



h(q\ 

miUi 

(k-q) 



(38) 



(39) 
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+(V S - v v 



"1 (h(k)-^ - 3/ B (A0l - miWi f/ 4 ^)^ + 3/ 6 (A;) 



M 2 

-|/3W^^-/ 5 (% 2 



M 2 



(40) 



Now the normalization condition for the wave functions with the components /3(g), /4(g), 
/5(g) and fe(q) is read as follows: 



dq I6U1UJ2 
(2tt) 3 3 



M 2 

3/ 5 / 6 - + 



2mia;i 2m\ijj\ 



2M. 



(41) 



Thus the results about the mass spectra and the wave functions for the J PC = 1 bound 
states can be obtained by solving the coupled Eqs.( l37H40|) numerically. 

II. The bound states with quantum numbers J PC = 2 ++ which in non-relativistic frame- 
work are 3 P2 and/or 3 F2 states mainly: 

The general form of the wave function for J PC = 2 ++ states can be written down as: 



AT 



M 2 



Mf 5 (q) + JPf 6 {q)] + 1 jh{q)f afh P<A±rm* 



(42) 



where the (symmetric in fi and v) is the tensor polarization of the meson. From the last 
two equations of Eq.Q: 

< 2 + + - + (g) = <- + + (g) = o, (43) 

we obtain the constraints on the components of the wave function: 



h(q) 



-q 2 h(q) + M 2 f 5 (q) 



h{q) 



U{q)M 



Mm\ mi 

Put the constraints into Eq. (l4"2"j) . then only four independent components /3(g) /4(g), /5(g) 
and fe(q) are left: 

_l q 2 



M Mmi 



)h{q) + ^tfM 



+(y + ^)Mh{q) + rffe(q) + i J -^e^P a q^ lll5 \ . 
nil mi J 



Mq) 



(44) 



From the formulation Eq. ff44|) . it is easy to see that the terms, which is proportional to /s or 
fe and with the factor (e^ u q l ['j^)M or (SuvQxl* 1 ) ? 1 are °f -P-wave nature (linear in gj_), and 
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the terms, which proportional to f'3 or f 4 and (Suv^J-IDm or (^vQxQD^p > are °f -F-wave 
nature (cubic in q±). Therefore, Eq. (j4*4"]) describes P — F wave mixing properly. 

Put Eq. fj44|) into the first two equations of Eq.Q and take various traces for 7-matrix 
on both sides of the equations, we obtain the coupled equations for the four independent 
components fs, f±, /s and f§ as follows: 



(M - a*) [fMj^ - fM) + [fMjf, + / 6 (g)) ^ 



d 3 k 



M 2 

1 ( -(Vs + V v ) q 2 U(k)^-f 5 (k)\ [k 2 q 2 -3(k-q)' 



(27r) 3 2o; 1 V ^ 

+mi(v; - V v ) 



. n .k 2 q 2 -?>{k-q) 2 , M 2 

-wi I /4W ^ - 2/ 6 (% 2 



(45) 



(M + 2a*) I (f 3 (q)£- 2 - /5(g)) - (fM^- 2 + /e(g)) ^} 

^ 1 [ {V s + K)g 2 (/s(*)^ - h(k)) [*V - 3(fc • q) 2 



(2tt) 3 2w 2 g 4 \ 
+mx(K-K) 



M 2 



-^1 ( Ai^j jf 2 2/ 6 (A;)g I 



(46) 



(M-2o; 1 ){-/ 5 (g)m 1 + / 6 ( g )a; 1 } 



d 3 £ 1 
(27r) 3 w 1 g 2 



x 



~2 (V. + K)/ 6 (*0 [A;V - 3(fc ■ g) 2 ] + (V s - V v )m 1 [(f 5 (k)u x - f 6 (k) mi ) 



- (/ S (*)wi + h{k) mi ) ^ + {h{k)u x + n{k) mi ) 



{k-q)\- 



(47) 



(M + 2wi) {-/ s (g)m! - / 6 (g)a;i} 



d 3 £ 1 
(2vr) 3 u;ig 2 



1 



x { ^(K + Ki)/ 6 (A;) [fc 2 g 2 - 3(A • g) 2 J - (K - K)™i [(/b^Wj + h(k) mi ) 



k 2 (k ■ q) 2 

Us{kpi - /i(fe)mi) + (/ 3 (fc)wi - hi^mx) 2 



(48) 
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TABLE I: Parameter of Vq in unit of MeV. 





cc 


bb 


nJ PC = nO- + ( 1 S ) 


-0.314 


-0.240 


nJ PC = n 1— ( 3 5i) 


-0.176 


-0.166 


nJ PC = n 0++( 3 P ) 


-0.282 


-0.174 


nJ PC = n l ++ ( 3 Pi) 


-0.162 


-0.141 


nJ PC = n 2++( 3 P 2 ) 


-0.110 


-0.121 


n J PC = n l + -( x Pi) 


-0.144 


-0.135 



Now the normalization condition is read as: 

/ wM: f™ 2 + 2hhf ~ 2,2/3 ( A £ + '•) } ^ 2M - < 49 > 

By solving the coupled equations Eqs. (j45N4"81) numerically, we obtain the mass spectra 
and relevant B.S. wave functions for the J PC = 2 ++ bound states. 

In fact, with the way described here, one may derive the BS-In equation for the other 
possible J PC states into their BS-CoEqs according one's wish. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section we solve the equations (BS-CoEqs) numerically and discuss the obtained 
results. 

Since the coupled integration equations are quite complicated, so we solve them numeri- 
cally only, and additionally with certain approximation such as that a cut on the up-bound 
of the integrations in the equations has been made. 

To solve the equations, we also need to fix the parameters appearing in the kernel Eq. (fl3|) 
although the kernel is based on QCD inspirer and the Cornell potential for non-relativistic 
heavy quark model as reference. Usually, the parameters are fixed by fitting the best exper- 
imental data. Since now quite a lot of data about the charmonium and bottomnium with 
quantum data J PC = 1 are available and quite precise, so the most parameters are fixed 
by the data. Since Vq in the kernel originates from QCD non-perturbative effects, its value 
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TABLE II: Mass spectra of (cc) and (bb) systems with quantum numbers J PC = + , 1 . Here 
( 2S + 1 )Lj denotes the dominant component in the state respectively. 'Ex' means the experimental 



results from PDG 



3 



(and the data for % come from reference 



n J PC (( 2S +VLj) 


Th(cc) 


Ex(cc) 


Th(66) 


Ex(66) 


1 O^So) 


2980.3(input) 


2980.3 


9390.2(input) 


9388.9 


2 (T+eSo) 


3576.4 


3637 


9950.0 




3 0-+(^ ) 


3948.8 




10311.4 




1 1-( 3 S!) 


3096.9(input) 


3096.916 


9460.5(input) 


9460.30 


2 1— ( 3 Si) 


3688.1 


3686.09 


10023.1 


10023.26 


3 1— ( 3 £>i) 


3778.9 


3772.92 


10129.5 




4 1"-( 3 5i) 


4056.8 


4039 


10368.9 


10355.2 


5 1— ( 3 £>0 


4110.7 


4153 


10434.7 




6 1— ( 3 5 x ) 


4329.4 


4421 


10635.8 


10579.4 


r I" ( 3 Si) 


4545.9 




10852.1 


10865 



is to account the states with various so we fix it by fitting the mass of the ground 

states. Thus the parameter Vq vary with J PC . 

By fitting data, the values of the parameters for all of the states are those as follows: 

a = e = 2.7183 , a = 0.06 GeV, A = 0.21 GeV 2 , 

and m c = 1.62 GeV, m b = 4.96 GeV . (50) 

Since the running coupling constant is used, so we also need to fix Aq C d- There are 
three active flavors for (cc) system, i.e. Nf = 3, accordingly we adopt Aqcd = 0.27 GeV 
and the coupling constant at the scale of charm quark mass, a s (m c ) = 0.38. There are four 
active flavors for (66) system, i.e. Nf = 4, so Aqcd = 0.20 GeV, and the coupling constant 
a s (mfc) = 0.23. By the fitting ground state data mainly, the fixed value of Vq for various 
J PC states is listed in TABLE I. 

The spectrum of charmonium and bottomnium (ground states and excited states) ob- 
tained by solving the coupled equations numerically is shown in TABLE II and TABLE 
III. Note here that since the couple-channel effects for the states above the threshold of 
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TABLE III: Mass spectra of (cc) and (bb) systems with quantum numbers J PC = ++ , 1 ++ , 2 ++ , 
l" 1 in unit of MeV. Here ( 2S+1 ^Lj denotes the dominant component in the state respectively. 
'Ex' means the experimental results from PDG 



n J PC ( 2S +^Lj 


Th(cc) 


Ex(cc) 


Th(66) 


Ex(bb) 


1 0++( 3 P ) 


3414.7(input) 


3414.75 


9859.0 


9859.44 


2 0++( 3 P ) 


3836.8 




10240.6 


10232.5 


3 0++( 3 P ) 


4140.1 




10524.7 




1 l++( 3 Pi) 


3510.3(input) 


3510.66 


9892.2 


9892.78 


2 l++( 3 Pi) 


3928.7 




10272.7 


10255.46 


3 l++( 3 Pi) 


4228.8 




10556.2 




1 2++( 3 P 2 ) 


3556.1(input) 


3556.20 


9914.4 


9912.21 


2 2++( 3 P 2 ) 


3972.4 




10293.6 


10268.65 


3 2++( 3 F 2 ) 


4037.9 




10374.4 




4 2++( 3 P 2 ) 


4271.0 




10561.5 




1 l+'CPi) 


3526.0(input) 


3525.93 


9900.2 




2 1+-{ 1 P 1 ) 


3943.0 




10280.4 




3 l+^Pi) 


4242.4 




10562.0 





'open-charm' or 'open-bottom' respectively have not been taken into account, so the results 
in TABLE II and TABLE III above the threshold cannot compare with experimental results 
directly. 

From the tables, one can read out the fine and hyperfine splitting precisely which are 
caused by the kernel with the fixed parameters. Therefore not only the gaps among the 
excited states and the ground states with fixed J PC but also the fine and hyperfine splitting 
among the states with different J PC are serious tests of the kernel and the B.S. approach to 
the heavy quarkonia. 

We cannot show all the numerical results of the wave functions which we have obtained 
here, alternatively, as typical examples, we only show some of the obtained wave functions 
with different quantum numbers J PC in figures FIGs.l,- • -,8. 

As usual cases, from the number of nodes of the wave functions in the figures we can 
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|q| GeV 



1.2 - 




|q| GeV 



FIG. 1: The solutions for the wave functions of the ground and the first excited state (from left to 
right) with quantum number J PC = h . The wave functions (solutions) of the low- lying states 
(the ground and the first excited state) with quantum number J PC = K The above two are 
those for charmonium and the below two are those for bottonium. 



realize how high an excited one or the ground each one of the obtained wave functions is. 

From FIG.l, we see clearly that for the states J PC = h the solution has the property 
fx — f%, so we may re- write the wave function Eq.( fl6|) as 



+ — )JPfi(q)+Mf 2 {q) 
mi 



'l + —)F + M 
mi 



7 5 



(51) 



here 4> Q _ + (q) — fi(q) — f2(q), the numerical solution of EqtflTl). 

From FIGs.2,3,4, we may see that the situations for the states with quantum numbers 
J PC = l + ~ J PC = ++ and J PC = 1 ++ are similar. Indeed the wave functions of the 
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|q| GeV 



|q| GeV 



25 -, 




FIG. 2: The solutions for the wave functions of the ground and the first excited state (from left to 
right) with quantum number J PC = l" 1 . The wave functions (solutions) of the low-lying states 
(the ground and the first excited state) with quantum number J = 1 H . The above two are 
those for charmonium and the below two are those for bottonium. 



ground state and the excited states for J PC = 1 + Eq.f l2~Tj) become 



= q± -e] 



1+ _(g)(g±-el) 



TOi ' M 



nix' M 



7s 



75 



(52) 



here <fi + _(q) — fi{q) — f2{q), the numerical solution of Eq. (|23l) ; the wave function for 
j pc = ++ Eq. (|25jl becomes 



(f p 0++ (q) = h{q)4± + h{q) 



4± F A(q)q 



7& 



0++ 



M 
M m 1 



rrii 



l± 



(53) 



19 





FIG. 3: The solutions for the wave functions of the ground and the first excited state (from left to 
right) with quantum number J PC = ++ . The wave functions (solutions) of the low-lying states 
(the ground and the first excited state) with quantum number J PC = ++ . The above two are 
those for charmonium and the below two are those for bottonium. 

here <p Q++ (q) — fi(q) — — f 2(0)1 the numerical solution; the wave functions for J PC = 1 ++ 
Eq.f l31|) become 



V? ++ (g) = ie^ aP P v qle%{q)M^ + f 2 (q)f>f + ih{q)e^ S q Lp P alslh /m 1 ]/M 



r2 



{q)e^ aP P u qle^ + + ie^ 5 q^ p P tjl5lb /m 1 /M 2 , (54) 



here (f> l++ (q) — fi(q) — —f2(q), the numerical solution of Eq.( |32|) . 

Moreover from FIG. 5 and FIG. 6 one may see the S — D wave mixing for the J PC = 1 
states, and from FIG. 7 and FIG. 8 one may see the P — F wave mixing for the J PC = 2 ++ 
states clearly. For the J PC = 1 states, from the figures (FIG. 5 and FIG. 6) we can see that 
for the first two states (the solutions for the ground one and the first excited one) the S- 
wave components /s and fa are dominant, and for the third state (the solution for the second 
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Ill GeV |q| GeV 

FIG. 4: The solutions for the wave functions of the ground and the first excited state (from left to 
right) with quantum number J PC = 1 ++ . The wave functions (solutions) of the low-lying states 
(the ground and the first excited state) with quantum number J PC = 1 ++ . The above two are 
those for charmonium and the below two are those for bottonium. 

excited one) the D-wave components f 3 and f± are dominant, etc. Therefore in TABLE I, 
we denote the first two states as J ( 3 Si) = 1 and the third one as J PC ( 3 Di) = 1 , 
etc in turn for high excited states. For the J PC = 2 ++ states, similarly we can see from 
the figures that the P-wave components f§ and /g are dominant in the first two states (the 
ground one and the first excited one), and the F-wave components f% and are dominant 
in the third state (the second excited state), etc. Therefore in TABLE II, we denote the 
first two states as J PC '( 3 P2) = 2 ++ and the third one as J PC ( 3 F2) = 2 ++ etc, in turn for 
high excited states. 

Of the J PC = 1™ states we may also see from the figure that for the .S-wave dominant 
states they have the properties: f 3 ~ — / 4 = (q) and f 5 ~ — / 6 = if) (q), so the 
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FIG. 5: The wave functions (solutions) of the low-lying states, the ground and the first four excited 
states, (from left to right) for charmonium with quantum number J PC = 1 . 




FIG. 6: The wave functions (solutions) of the low-lying states, the ground and the first two excited 
states, (from left to right) for bottonium with quantum number J = 1 . 

solutions (wave functions) can be re- written (from Eq.( l36|) ) as 



Mrm M ' M 2 
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FIG. 7: The wave functions (solutions) of the five low- lying states, the ground and the first four 
excited states, (from left to right) for charmonium with quantum number J PC = 2 ++ . 




FIG. 8: The wave functions (solutions) of the five low-lying states, the ground and the first four 
excited states, (from left to right) for bottonium with quantum number J PC = 2 ++ . 

+ t- fa) Wx + *x • <>A) - [AF + - ^ } ; (55) 

[ V 777]/ L 777,1 777! J J 

whereas of the .D-wave dominant states they have / 3 ~ / 4 = (g) and f$ ~ — f$ = 
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■0 (q), so the solutions (wave functions) can be re- written as 

,2 



± 



<i>i—(q)(qi 



M> M 2 



M 



(56) 



+ (q)\ (m£ + q ± ■ el—) + 

From the figures (FIG. 7 and FIG. 8) of the J PC = 2 ++ states, we may see that for 
the P-wave dominant states, the solutions have the properties: /s ~ —fa = -0 2++ (q) and 
h - -h 



2 ++(q) an d the solutions (wave functions) can be re-written (from Eq.(jHJ) as 

W Mrm 



2+" 



Xq>*q±qt 



-H0. 



2+- 



p I mi mi 



M 2 



(57) 



whereas for the F-wave dominant states, the solutions have the properties: / 5 ~ /@ = 
■0 ++ (g) and / 3 ~ / 4 = 2++ (g) and the solutions (wave functions) can be re- written (from 
Eq.flBD) as 

\.d. a 2 

(q±) ^ <f> 2++ (q)e*qlq± 



P,2++ 



(ILL 



Mmi 



+ 



M 2 



2 + 4 



(58) 



^1 ITll ) 

Finally we would like to discuss the wave mixture further. As shown in TABLE. II and 
Eqs. (l55|l56l) . each of the states (either charmonium or bottomonium) for J PC = 1 contains 
S — D wave mixing. Some are S'-wave dominant and the others are P-wave dominant. The 
third one of charmonium with m = 3778. 9MeV below the threshold of 'open-charm', for 
instance, clearly is a P-wave dominant state, and in its decay into II, (I = e,fi) only its 
.S'-wave components play a role, so the fraction width of the pure leptonic decay is compara- 
tively small. Indeed it corresponds to the observed one ■0" (with mass m = 3772.92MeV) as 
pointed in 2j. Whereas, similarly there is wave mixture for bottomonium too, for instance, 
once more the third one of bottomonium with m = 10129. 5MeV below the threshold of 
'open-bottom' is also a P-wave dominant state and it decays into II, (I = e,/i) only via 
its S'-wave components, so the fractional width of the pure leptonic decay is comparatively 
small too 3 . Furthermore, the fractional width of the pure leptonic decay for such a state will 



3 In Rcf . [2] , the wave mixture is obtained via additional interaction and that of charmonium is concerned 
only, i.e. the mixture for bottomonium is not discussed, that is different from here. Here the mixture for 
charmonium and bottomonium is fully determined by the B.S. kernel well, so we need to discuss the sense 
for bottomonium on experimental observations too. 
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be comparatively much smaller (a quarter) than that of charmonium, due to the two factors 
that the charge of bottom-quark is smaller than that of charm-quark, and the compara- 
tive weight of the S'-wave component to the D-wave component in the D-w&ve dominant 
state, which is proportional to v 2 {v bot tomonium < v charmonium ) as indicated in Eq.([36])), is 
small. Therefore such a state is very difficult to be observed either in e + e~ energy scanning 
experiments at CLEO and B-factories (due to low production rate) or in hadron colliders 
(due to very small branching ratio for the lepton pair decay and various backgrounds etc). 
We believe that all such .D-wave dominant states for bottomonium must be still missing 
in experiments so far, even if our prediction here is true. We conjecture that such states 
may be observed at Z-factory such as Giga-Z etc elsewhere via e + e~ — > {bb)\ — (3 Dl ) + 7 or 
e + e~ — > (bb)\ — (3D!) + • • •, because there the backgrounds can be controlled comparatively 
easy, and numerous such bottomonium states enough for experimental observation can be 
produced via on-shell Z-boson [20]. For the P — F wave mixture, since the first state with 
quantum numbers J PC = 2 ++ is a 'high' excited state already so there are only fewer of 
the J PC = 2 ++ states below the threshold of 'open-charm' or 'open-bottom', thus there are 
fewer example states which can be used to test the wave mixture, although the tests and 
the situation essentially are quite similar to the cases of 1 states for S — D wave mixture. 

When the bound states does not consist of a pair of quark and antiquark (not as charmo- 
nium and bottomonium here), the quantum number C is not a good one, then the present 
way to solve the problem (the relevant B.S. equation) should be changed accordingly, but 
its main steps may be still kept and interesting results, which are different from the present, 
are obtained finally. In fact we have considered the double heavy system (cb) or (cb) as an 
example for non-(quark-antiquark) binding system and solved the relevant B.S. equation in 
a similar way, but due to differences we put the results and discussions about the double 
heavy system (cb) or (cb) elsewhere in Ref. 211 ]. 
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